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ABSTRACT 

The central parsec of the Galaxy contains dozens of massive stars with a cumulative mass loss rate 
~ 10~^ Mq yr~^. Shocks among these stellar winds produce the hot plasma that pervades the central 
part of the galaxy. We argue that these stellar wind shocks also efficiently accelerate electrons and protons 
to relativistic energies. The relativistic electrons inverse Compton scatter the ambient ultraviolet and 
far infrared radiation field, producing high energy 7-rays with a roughly constant luminosity from ~ 
GeV to ~ 10 TeV. This can account for the TeV source seen by HESS in the Galactic Center. Our 
model predicts a GLAST counterpart to the HESS source with a luminosity of « 10'^^ ergs s"^ and 
cooling break at w 4 GeV. Synchrotron radiation from the same relativistic electrons should produce 
detectable emission at lower energies, with a surface brightness w 10'^^i?^3 ergs s~^ arcsec"^ from ~ 
THz to ^ keV, where i3_3 is the magnetic field strength in units of mG. The observed level of diffuse 
thermal X-ray emission in the central parsec requires B <, 300 fi G m our models. Future detection of the 
diffuse synchrotron background in the central parsec can directly constrain the magnetic field strength, 
providing an important boundary condition for models of accretion onto Sgr A*. 

radiation mechanisms: thermal - radiation mechanisms: 



Subject headings: Galaxy: center 
non-thermal 

1. INTRODUCTION 

The young stellar cluster in the central parsec of the 
Galaxy contains nearly two dozen luminous blue super- 
giants and Wolf-Rayet stars, in addition to a larger pop- 
ulation of lower mass O and B stars (Krabbe et al. 1995; 
Genzel et al. 2003). These stars were produced in a short 
lived compact starburst that occurred ~ 10^ — 10^ years 
ago. Collectively, the massive stars in the central par- 
sec lose w 10~^ M0 yr^^ in stellar winds with velocities 
ranging from ~ 300 to 1000 km s"^ (Najarro et al. 1997; 
Paumard et al. 2001). Shocks produced by colliding stel- 
lar winds can account for the hot (ss 1 — 2 keV) X-ray 
emitting plasma that fills the central few parsecs of the 
Galaxy (Quataert 2004; Rockefeller et al. 2004; Cuadra 
et al. 2005). These stellar winds are also the most likely 
source of fuel for the central massive black hole Sgr A*: 
~ 0.3 — 1% of the mass lost to stellar winds is gravita- 
tionally captured by the black hole and flows in to smaller 
radii (e.g., Coker & Meha 1997; Quataert 2004; Cuadra et 
al. 2005). The remaining mass is thermally driven out of 
the central stellar cluster in a wind. 

In this Letter we investigate the possibility that stellar 
wind shocks in the Galactic center also efficiently accel- 
erate electrons and protons to ultra-relativistic energies. 
This claim is motivated by the analogy between the strong 
collisionless shocks expected in the Galactic center and 
similar shocks in supernovae remnants and 7-ray bursts 
that are inferred to be efficient sites of particle acceler- 
ation (e.g., Drury 1983; Blandford & Eichler 1987). We 
show that a population of relativistic particles in the cen- 
tral parsec of the galaxy produces detectable emission from 
the radio to the 7-rays (§2). We then compare the predic- 
tions of our model with current observations and discuss its 



further implications (§3). We focus primarily on emission 
by relativistic electrons because it has fewer uncertainties, 
but we also briefly consider pion production by relativistic 
protons. 

2. PREDICTED SPECTRUM 

The total kinetic energy thermalized by shocks in the 



Galactic center is E ^ 
-1 ,,tU«,.« A/f — in-3 
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« 10^«/o.3M_3«8^ ergs 

s"", where = IQ^'^M-^ Mq yr~^ is the total stellar 
mass loss rate, = lOOOuskms"^ is the average wind 
speed, and / = O.3/0.3 accounts for geometric uncertain- 
ties such as the fraction of the stellar wind material that 
undergoes strong shocks (the normalization for / is taken 
from the 3D simulations of Rockefeller et al. 2004) . 

We assume that a fraction 77 = 10~^7y_2 of the shock en- 
ergy goes into accelerating relativistic electrons, where the 
normalization of rj is motivated by observations of super- 
nova remnants (see, e.g., §2.2.2 of Keshet et al. 2003 for 
a summary of the relevant observations) . Thus the energy 
supplied to relativistic electrons is 

Ee « 1036/o.3?7-2A/-3«8 ergss-^. (1) 
We assume for simplicity that the electrons are injected 
with a spectrum oc 7"^ with p = 2, similar to the index 
typically observed in supernova remnants (e.g., Brogan et 
al. 2005; Aharonian et al. 2005) and theoretically ex- 
pected for strong shocks (Bell 1978; Blandford & Ostriker 
1978). 

The high energy electrons cool by synchrotron radiation 
and inverse Compton (IC) emission. Figure 1 summarizes 
our predictions for the spectrum produced by these two 
processes; the main features in this spectrum are described 
analytically below. The magnetic field in the central par- 
sec is poorly constrained. We normalize our estimates us- 
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ing B = i?_3 mG, which is cquipartition with the hot ther- 
mal plasma and is similar to the field strength inferred us- 
ing the properties of non-thermal radio filaments on larger 
scales in the Galactic center (e.g., Morris & Serabyn 1996; 
see, however, LaRosa ct al. 2005). Equipartition field 
strength could in principle also be generated by the shocks 
themselves (e.g., Medvedev & Loeb 1999). 

The young stellar cluster in the central parsec has an ul- 
traviolet (UV) luminosity of L = 10^^L4i ergs s~^ within 
a radius R = 0.5-Ro.5 P*^ (Krabbe et al. 1995), implying 
a radiation energy density of Uph = lO"''?/^^ ergs cm~^ 
with a characteristic photon energy of Eph ~ 10 eV. Only 
w 20% of the UV radiation is absorbed by dust in situ 
(Davidson et al. 1992), which leads to a far infrared (FIR) 
radiation field with £ph ~ 0.04 eV and U}^ ~ 0.1 or equiv- 
alently i/i^ ~ 1. 

The relativistic particles in the Galactic center are ad- 
vected out of the central parsec with the thermally driven 
wind on a timescale texp ~ R/vw « 500 R0.5/VS years. where 
Sufficiently relativistic electrons with 7 > 7c cool on a 
timescale <^ t^xp- For our fiducial normalization the cool- 
ing is primarily by IC radiation (unless B ^ 1.6 mG). In 
this case the critical Lorentz factor where the IC cooling 
time equals the expansion time is given by 

^^^STrmeC^^^xlO^^^ (2) 
(TtL L41 

where ctt is the Thompson cross section. Electrons with 
7 > 7c lose all of their energy to IC radiation before es- 
caping the central parsec, producing a power 



6 X 103Vo.3??-2M-3f| ergs s'^ (3) 



where 7^02; is the maximum Lorentz factor of the acceler- 
ated electrons. The exact value of ^max is uncertain but 
a rough estimate can be obtained by balancing the cool- 
ing time with the acceleration time at the shock front, 
tacc ~ ^~^(^ hwi where f2 is the relativistic cyclotron fre- 
quency (e.g., Loeb & Waxman 2000). This yields 

1 1/2 



eB 



— ~3xlO^^LJ. (4) 
c 



Note that in the above numerical estimate we have nor- 
malized the photon energy density using the FIR radiation 
field, while in our estimate of 7c in equation (2) we used 
the full UV radiation field. This is because electrons with 
7 ^ "/max see UV photons in their rest-frame as 7-rays 
with energy ^ m^cP and do not {ffic;iently scatter the UV 
photons due to the Klein-Nishina suppression of the cross 
section. 

Equations (2) and (4) imply that the IC power in equa- 
tion (3) applies for photon energies Emin < E < Emax, 
where 

^ ^ Eph 



Err 



%Eph ~ 4 GeV 



and 



Emax « Imax^ph « 30 TcV 
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(5) 



(6) 



Relativistic corrections to the electron scattering cross sec- 
tion lead to an additional suppression of the IC power 
above 

Eph 



Eb « 25 TeV 



0.04 eV 



(7) 



For E < Emim the electrons do not lose all of their en- 
ergy before flowing out of the central parsec and so the IC 
power is suppressed, with vL^ oc {E / EminY^"^ ■ 

2.1. Synchrotron Emission 

The same electrons that produce IC power in the 7-rays 
will also produce synchrotron emission at lower energies, 
with Pa/ Pic = Us/Uph- For i/^ < 1^ < i^max the total 
synchrotron power is (assuming Ub < Uph) 

52 
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fo.3V-2M.3vl ergs s 1 (8) 
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(9) 



(10) 



Near Vmax the electrons only IC scatter FIR photons and 
so Uph ~ should be used in equation (8).^ 

Equations (8)-(10) imply that synchrotron emission 
from shock accelerated electrons produces a diffuse radia- 
tion field from the radio to the X-rays in the central par- 
sec. This emission would be extended on a radial scale of 
» 0.5i?o 5 pc i=a 10" and so would have a surface brightness 
of 

7, » 7 X 10^^ I I fo.3V-2M_3vl ergss"^ arcsec'^ 



L41 



2.2. -y-ray Emission from Pion Decay 



(11) 



Relativistic protons can also be a significant source of 
high energy emission through the production of neutral 
and charged pions; here we briefly estimate the expected 
7-ray luminosity from neutral pion production in the cen- 
tral few parsecs of the Galaxy. We assume that a larger 
fraction, ^ 10%, of the shock energy goes into accelerating 
relativistic protons (Blandford & Eichler 1987), in which 
case 

(12) 



Ep ~ 10-"fo.3ii-iM_3vieTgss- 



The hot thermal plasma in the Galactic center has a 
density of n w 30 cm~^ (e.g., Baganoff et al. 2003). 
The timescale for relativistic protons to lose their en- 
ergy via neutral pion production is tpion ~ 10^ years. 
Comparing this to the expansion time of the hot plasma 
~ 500 years, we find a very low GeV 7-ray luminos- 
ity from pion production on the ambient hot plasma, 
uL^ « 10-'^Ep/ln{^rnax) ~ 1032ergss-i. 

A more promising source of 7-rays from pion decay in 
the central few parsecs is the interaction between the shock 
accelerated protons and the circum- nuclear disk (CND). 



U}^s J V0-04eF, 

^ Note that if ~ 1 and B_3 ~ 1, Ub > C/ph and so equation (8) would not apply. As we discuss in the next section, however, observations 
constrain B ^ 300//G so this possibility is not realized. 
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The CND is a ring of dense molecular gas located « 1.5 
pc from Sgr A*, with a scale height of ~ 0.5 pc (Jackson 
et al. 1993; Christopher et al. 2005). The mean density of 
the CND is at least n « lO^ns cm~^, and Christopher et al. 
(2005) argue that it could be several orders of magnitude 
higher. The morphology of the diffuse X-ray emission in 
the central parsec suggests that the outflowing wind from 
the central star cluster impinges on and is confined by the 
CND (Baganoff et al. 2003; Rockefeller et al. 2004). 

To estimate the 7-ray luminosity from pion decay in the 
CND, we assume that a fraction ^ = O.l^o.i of the rel- 
ativistic protons created in the central 0.5 pc enter into 
the CND, and that a fraction ^ 1/3 of the proton energy 
goes into 7-rays, rather than neutrinos or other particles. 
The normalization for ^ is based on the solid angle sub- 
tended by the inner edge of the CND, but in principle 
^ could be much smaller if relativistic particles are mag- 
netically shielded out of the CND. The cooling time for 
relativistic protons due to pion production in the CND is 
tpion ~ 10^ n^^ years. By contrast, the timescale for rela- 
tivistic protons to escape from the CND is uncertain, and 
depends sensitively on the magnetic field geometry and 
cross-field diffusion. Unless the effective diffusion velocity 
is ^ 1000 km s~^, however, the protons will lose all of their 
energy in the CND because the cooling time is shorter the 
escape time. Under this assumption, the 7-ray luminosity 
from the CND is 

i/L, - ~ 2 X 103Vo.3eo.ir?-iM-3«|ergss-i. 

Equation (13) predicts the maximum 7-ray luminosity not 
far below that produced by relativistic electrons (eq. [3]), 
but the estimate is much less certain because the relativis- 
tic protons might quickly leak out of the CND along open 
magnetic field lines. The spectrum of 7-rays from pion 
decay would extend from w 70 MeV to many TeV, with 
the maximum photon energy determined by the maximum 
energy of the accelerated protons. 

3. DISCUSSION 

Figure 1 summarizes our estimates of the IC and syn- 
chrotron radiation expected from shock accelerated elec- 
trons in the central parsec. We do not include the pre- 
dicted radiation from pion decay estimated in §2.3 because 
it is much less certain, though we reiterate that the 7-ray 
flux from pion decay in the CND could be non-negligible. 

Figure 1 also summarizes some of the data relevant to 
our predictions. In particular, IC radiation of the ambi- 
ent FIR radiation field can account for the TeV source 
seen by HESS (Aharonian et al. 2005) and CANGAROO 
(Tsuchiya et al. 2004) towards the Galactic center. In 
Figure 1 we have adjusted the normalization of the IC 
component to account for the HESS luminosity by setting 
?7 « 5 X 10~^ in equation (3). Our model naturally ac- 
counts for the HESS spectrum and makes the strong pre- 
diction that the Gamma-Ray Large Area Space Telescope 
{GLAST)'^ should observe a cooling break at ~ 4 GeV. 
The energy of the cooling break is very well determined in 
our model (eq. [5]) since it is set by the UV radiation en- 

* See http://glast.gsfc.nasa.gov 

^ On much larger (> 2') scales, LaRosa et al. (2005) argue that the fi 
has been used to argue for B ~ mG (e.g., Morris & Serabyn 1996). 



ergy density and the expansion time of the thermal plasma, 
both of which are observationally well constrained. 

Atoyan & Dermer (2004) also explain the TeV emission 
from the Galactic center as IC scattering of FIR photons 
by relativistic electrons (sec Pohl 1997 and Mclia et al. 
1998 for related ideas in the context of the EGRET source 
3EG J1746-2851, which is now beUeved to be well-offset 
from the central parsec; Pohl 2005). In their model, how- 
ever, the electrons are produced in the termination shock 
of an outflow from Sgr A*. On energetic grounds, the 
shocked stellar winds investigated in this Letter may be a 
more important source of particle acceleration: the accre- 
tion rate onto Sgr A* is believed to be ^ 1O^*M0 yr"-*^ 
(e.g., Agol 2000; Quataert & Gruzinov 2000), implying a 
total accretion power « O.lMc^ ^ 5 x 10^'^M„8 ergs s~\ 
somewhat smaller than the energy thermalized by shocked 
stellar winds in the central parsec. In addition, the pres- 
ence of stellar wind shocks in the central parsec is ob- 
servationally well established, while the properties of an 
outflow from Sgr A* (its mass loss rate, collimation, etc.) 
are significantly less well understood. 

In addition to IC emission at high energies, we also 
predict that synchrotron radiation from shock accelerated 
electrons in the central parsec leads to a diffuse back- 
ground of radiation from ~ THz to ~ keV (Fig. 1), 
with a surface brightness given by equation (11). This 
emission would be extended on the scale of the central 
stellar cluster, « 10". The magnitude of the diffuse syn- 
chrotron emission is sensitive to the magnetic field strength 
in the central parsec, which is poorly constrained ob- 
servationally'' The most stringent observationally limits 
to date on diffuse emission at the relevant wavelengths 
are probably from deep Chandra observations which find 
extended thermal emission with a surface brightness of 
« 7.6 X 10^^ ergss"^ arcsec"^ in the central few parsecs 
(Baganoff et al. 2003). Using equation (11) with L41 « 0.1 
appropriate for the FIR radiation field scattered by elec- 
trons with 7 ~ 7max, we find that the Chandra observa- 
tions require that the magnetic field strength in the central 
parsec be ^ 300/iG in order to not overproduce the ex- 
tended X-ray emission. It is worth highlighting that future 
probes of the nonthermal diffuse emission could directly 
constrain the magnetic field in the central parsec. In par- 
ticular, the ratio of the luminosity of the diffuse emission 
to the luminosity of the HESS TeV source is determined 
by Ub/U^^^ and thus can be used to measure Ub- Such a 
measurement would provide a very useful outer boundary 
condition for models of accretion onto Sgr A*, which are 
sensitive to the geometry and magnitude of the magnetic 
field at large radii (e.g., Igumcnshchev & Narayan 2002). 

Lastly, wc note that some fraction of the shock acceler- 
ated electrons and protons in the central parsec should be 
gravitationally captured by Sgr A* along with the thermal 
plasma, providing a seed population of relativistic particles 
to the accretion flow. The Lorentz factor of such particles 
will, at a minimum, increase via adiabatic compression as 
7 oc p^/'^ . The density contrast between the Bondi accre- 
tion radius {Rb ~ IQ^Rs) and the event horizon is beheved 
to be a factor of ~ 10^ (see Quataert 2003 for a review), 

strength is ~ 10/iG, while the prevalence of nonthermal radio filaments 
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implying that mildly relativistic electrons at i?^ have 
7 ~ 100 close to the black hole. Such electrons would emit 
in the mm-MIR and could contribute significantly to the 
observed emission from Sgr A*. This possibility will be 
explored in more detail in a future paper. 
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Fig. 1. — Predicted emission from shock accelerated electrons in the central parsec of the Galactic center. The solid line shows the IC 
radiation (eq. [3]) while the dashed line shows the synchrotron radiation for B = 0.3 niG and Upi^ = 10^^ erg cm^'^ (eq. [8]). The parameter 
r) measuring the fraction of energy in shock accelerated electrons is adjusted to reproduce the HESS TeV luminosity (r; ft; 5 X 10^'^). With 
this normalization the 7-ray spectrum is uniquely determined, including the predicted cooling break in the GLAST band at ft; 4 GeV; the 
flux of synchrotron radiation depends primarily on the magnetic field strength in the central parsec. The high energy cutoffs to the IC and 
synchrotron radiation arc determined by "fmax, the maximum Lorentz factor of accelerated electrons (eq. [4]). Klcin-Nishina suppression of 
the electron scattering cross section produces a high energy cutoff to the IC emission at comparable energies (cq. [7]). The exact spectral 
shape above the high energy cutoff is uncertain and depends on the details of of the electron distribution around -ymax ■ The two power-laws 
shown for the HESS data indicate the uncertaintes in the power- law slope (Aharonian et al. 2004). 



